Mon. Not. R. Astron. Soc. 000. ITHTT1 (2009) Printed 13 September 2010 (MN LflftX style file v2.2) 



First experimental demonstration of temporal 
hypertelescope operation with a laboratory prototype 
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ABSTRACT 

In this paper, we report the first experimental demonstration of a Temporal Hyper- 
Telescope (THT). Our breadboard including 8 telescopes is firstly tested in a manual 
cophasing configuration on a ID object. The Point Spread Function (PSF) is measured 
and exhibits a dynamics in the range of 300. A quantitative analysis of the potential 
biases demonstrates that this limitation is related to the residual phase fluctuation 
on each interferometric arm. Secondly, an unbalanced binary star is imaged demon- 
strating the imaging capability of THT. In addition, 2D PSF is recorded even if the 
telescope array is not optimized for this purpose. 
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1 INTRODUCTION AND CONTEXT 
1.1 Hypertelescope concept 

High resolution Optical imaging instruments based on aper- 
ture synthesis have been developed over the last decades 
with the aim of reaching angular resoluti on in the nan o ra- 
dian range. These different instruments (|Lawsonlll9"9m use 
the property of the Zernike van Cittert theorem to recover 
the intensity distribution of the object by means of spa- 
tial coherence analysis. With this method, the instrument 
can never select the light coming only from one of the pix- 
els composing the full object because the measuremnts are 
being carried out on the Fourier spectral domain. For high- 
dynamics objects, such as a star + exoplanet system, this 
technique is limited since the information on the faint ob- 
ject is always mixed with the light emitted by the main 
source. Consequently, direct imaging is to be preferred and 
the analysis of the object is made easier in the image do- 
main than in the Fourier spectrum one. Since the beginning 
of high resolution imaging, measurements have never been 
achieved both with a very high resolution in the range of 
nanoradian and a very high dynamics in the range of 10 6 . 
In order to meet this challenge A. Labeyrie has proposed 
a sol ution which is known as the hypertelescope (|Labevrid 
1996). This new type of instrument solves the problem of the 
highly structured Point Spread Function (PSF) of a diluted 
array thanks to a pupil densification process. The PSF of a 
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hypertelescope being sharp and smooth, it is possible to use 
the instrument for direct imaging. The image / which equals 
the convolution of the object O by the PSF, looks like the 
object but with a limited resolution. Different versions of 
hypertelescopes ha ve been proposed using field combination 
in the pupil plane ([Vakili et all [2004) or pupil d ensification 
than ks to the use of monomode optical fibres ([Patru et al 
l20Q8h . 



1.2 Principle of classical Hypertelescope 

Parallel to the hypertelescope studies promoted by A. 
Labeyrie, we have proposed a temporal alternative to the 
initial design that used spatial classical optics (cf Fig[TJ). 
The main purpose of this new concept is to answer some 
technical difficulties met with classical hypertelecopes and 
to propose new functi o naliti es. 

In lRevnaud &; Delagd (|2007h we theoretically demonstrated 
the possibility to design a hypertelescope by using temporal 
optical path modulation. In the next paragraphs, we briefly 
recall the principle of a classical and a temporal hypertele- 
scope. 

Figure [T] recalls the structure of a hypertelescope as pro- 
posed by A. Labeyrie. This simplified drawing does not de- 
tail the reconfiguration and densification process. This tech- 
nique makes it necessary to remap the input pupil taking 
care to apply an homothetic contraction of beams center 
distribution in the output pupil according to the "golden 
rule of imaging interferometry" . Denoting 1/G the pupil re- 
configuration ratio, the distance between telescope Ti and 
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Figure 1. Spatial hypertelescope : the input pupil is reconfigured 
with a G contraction ratio and a densiflcation of the output pupil 
(TTi = G • OOi). The image is displayed as a function of spatial 
coordinates (x;y) in the image plane. 



the input pupil center T is related to the position OOi of 
the sub pupil i in the output pupil by: 

TTj 
G 

where O is the output pupil origin. After passing through 
the focusing lens, the beam i in the image plane generates 
a limited plane wave: 



OOi = 



(1) 



Ei- 



Ki 



• exp (ji/ji) • d(X) • exp (jKi • X) 



with 

_ 2tt OOj 



2tt TTj 
T ~Gf 



(2) 



(3) 



where / is the focal length of the focusing lens L, en 
and ^i are respectively the modulus and input phase of 
beam i. X is the position vector in the image plane with 
two coordinates x and y (X = (x;y) ). In the Labeyrie 
configuration d(X) denotes the field envelop resulting from 
the diffraction of each sub pupil supposed to be identical for 
all beams. The last term exp (jKi • X) describes the linear 
variation of the phase in the image plane with ol% and Pi 
slope coefficients along x and y axes. 



Ki = 27T(a z ,/3i) 



with 



(OOi) 
A/ 



and 



Pi 



(OOi) 
A/ 



(4) 



(5) 



where (OOi) x and (OOi) y represents the projection of vector 
along x and y axis. 

This way, each position TTi of telescope drives the slope 
(ai, Pi) of the corresponding wave front and reaches the im- 
age plane. The phase evolution ifi(x, y) = KiX can be anal- 
ysed along different axes parallel to the x axis. The phase 
variation ipi(x) is linear with an offset depending on the y 
coordinate : 

ipi(x) = 2tt • (ai • x + Pi • y) = 2tt • m • x + <Pio(y) (6) 
with 

(fio(y) = (fio(0; y) = 2ir- Pi-y (7) 

As demonstrated by A. Labeyrie, the possibility to re- 
trieve the image results from the coherent superposition of 



these different limited plane waves and the convolution of 
the resulting Point Spread Function (PSF) with the object 
intensity distribution. 

The PSF is given by the intensity distribution corre- 
sponding to the coherent superposition of field E{ in the 
image plane when the input pupil is illuminated by a plane 
wave with a wave vector perpendicular to the input pupil 
plane. Under such conditions, the ai coefficients are identi- 
cal (set to one for example) and the ^ phase terms in the 
equation [2] are equal to 0. 



PSF(X) 



d(X) 



^exp (jKi-X) 



(8) 



If in the Labeyrie configuration, d(X) corresponds 
to the sub pupil far field, in the IRAN (Interferometric 
Remapped Array Nulling) proposal this term is replaced by 
the sub-pupil function itself. In both cases this envelope term 
fixes the phase modulation span in the recombining plane. 
This phase range is directly related to the optical pieces di- 
mension and will not be easily adjusted in a real experiment. 

For a tilted point like source illuminating the telescope 
array with a — (9 X , y ) obliquity, the amplitudes ai remain 
constant but the phase ipi (as defined in eq J2J)) becomes: 



A 



(9) 



The total phase of the Ei field can be written in the 
image plane as : 



A A LrJ 



(10) 



This additional ipi term induces a spatial shift of the 
corresponding intensity distribution : 



h{X) =PSF(X + Gf0) 



(11) 



For an extended object, the incoherent superposition of 
the different contribution of the object O(0) leads to the 
fallowing image intensity distribution : 



/ 



I(X) oc / O{0) ■ PSF(X + Gf0)d0 



(12) 



Of course, any process able to achieve such phase modu- 
lations will provide an image with the same basic properties. 
The next paragraph deals with the possibility to obtain such 
a result in the time domain. 



1.3 From spatial to temporal domain 

As above mentioned, in a spatial hypertelescopes, the linear 
phase modulations Lpi(x) are related to the ai and pi tilts 
of the optical fields Ei reaching the observation plane. 

In the temporal case, the phase modulation <fi(t) is gen- 
erated thanks to an optical path variation induced by optical 
path modulators linearly actuated as a function of time as 
shown on figure [2] 



ipi(t) = 2tt • Vi • t + <Pio(y) 



(13) 



Using a temporal frequency Vi scaled on the a% slope 
allows to propose a temporal configuration equivalent to the 
spatial concept. For this purpose, the frequency Vi has to be 
proportional to ai with an arbitrary A— — coefficient. 
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Figure 2. Phase variation in the spatial and temporal cases : 
In the spatial case (top), the phase variation is induced by the 
obliquity of plane wave related to each aperture. In this case, the 
phase span is limited by the diffraction spot size. In the temporal 
case (bottom), a phase modulation is applied by the operator to 
temporally reproduced the slope observed in spatial case. The 
phase span is now driven by the optical path stroke and can be 
adjusted without any limitation. 
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According to this spatial to temporal transposition, the 
optical field: 

Ei = en - exp (jipi) • d(X) • exp {j[2n • on • x + y>io(y)]} (15) 

is replaced by : 

Ei = ai- exp (jipi) • d(t) • exp {j[27i • Vi • t + <Pio(y)]} (16) 

d(t) function expresses the span of phase variation. In the 
classical concept, the image is analyzed along the x axis. 
The A factor represents the relationship between the spatial 
parameter x and time t. 



x = A-t 



(17) 



In the time domain, the image is obtained by using the 
temporal PSF given by: 

2 



PSF(t;y) 



d(t) ■ exp {i[27r ■!/»•*+ ipio(y)]} 



(18) 



As seen in the previous case, for a tilted point like source 
with a = (0 X ; y ) obliquity, the total phase of the Ei field 
becomes: 



(19) 



The corresponding intensity distribution is a shifted tempo- 
ral PSF: 



Ie(At;y) = PSF[(At;y) + Gf9] 



(20) 



and the incoherent superposition of the different contribu- 
tion of the object O leads to the temporal image: 
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Figure 3. Temporal versus spatial display. The two-dimensional 
spatial display can be replaced by a temporal one. As long as 
time is a one dimensional parameter a set of scan is necessary to 
retrieve a time multiplexed image like in a video monitor. 



I (At; y)oc O(0) • PSF[(At; y) + Gf0]d0 



(21) 



This result demonstrates the full equivalence between 
spatial and temporal display as illustrated in figure [3] 

In order to analyze a two dimensional image, the (fio(y) 
phase has to be modified step by step in order to display suc- 
cessively the different s of the image using a set of parallel 
scans. The information provided along the x axis at a given 
position y is fully equivalent to the one observed as a func- 
tion of time t for the scan related to the same (fio(y) value. 
The image is temporally multiplexed as for a video raster 
scan. It can be noticed that in the spatial domain, the x span 
is determined by the diffraction field distribution related to 
the sub pupil geometry. In the temporal domain, the image 
field is determined by the extension of phase modulation di- 
rectly driven by the span of the optical path modulator. It 
can be adjusted very easily and reduced down to contrary 
to the spatial configuration in which the limitation results 
from the beam dimensions. The different advantages of tem- 
poral hypertelescopes are discussed in iRevnaud &; Delagel 
(2007|). 



2 IMPLEMENTING A TEMPORAL 
HYPERTELESCOPE 

2.1 General layout 

Figure H] shows a generic representation of a temporal hyper- 
telescope. This sketch uses optical fibres and an optical cou- 
pler but classical mirror trains, air delay lines and classical 
beam splitters could be used to design such an instrument. 
Let's consider a telescope array pointing at a scientific tar- 
get. The light is picked up at the N telescope foci and 
passes through optical path modulators before being recom- 
bined using an N to one beam combining coupler in a coax- 
ial configuration. The required optical path modulations of 
the N interferometric arms are generated using optical fibre 
stretchers. They induce an optical path variation Si in order 
to generate the convenient phase modulations as previously 
mentioned. For this purpose, the optical path modulation 
can be expressed as : 



A(OQ % 



(OOi) y • y 



(22) 



/ / 

These optical path modulations can be servo controlled 
using opto-electronic systems previously developed for such 
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Figure 4. Proposal for a temporal hypertelescope using optical 
fibres and a guided coupler acting as a combiner. The optical 
path modulators induce the phase variations necessary to form 
an image. The image is temporally multiplexed. 

kinds of applications (| Pelage et all (|2000h ; I Olivier et all 
(2005)). It allows to monitor, with a nanometric accuracy, 
the linearity of the optical path variation as a function of 
time. For each scan, the offset ^°°^ v y is set in order to 
display the signal that would be observed at y position for 
a classical spatial hypertelescope. 

To properly operate a temporal hypertelescope, the op- 
tical path modulators are driven by an 8 channel function 
generator and related high voltage electronics (not drawn 
in the picture). The output voltage drives the optical path 
modulators with a full span in the range of tens jim and a 
typical nanometric sensitivity. The electronic gain and the 
voltage generator slopes allow to set the Vi frequencies at 
the proper values. In such a configuration, we can theoreti- 
cally get the same imaging properties as for the first classical 
design using spatial pupil densification. The breadboard de- 
scribed in this paper and the related experimental results 
reported in the next paragraphs aim to demonstrate the va- 
lidity of this new concept. 




Figure 6. Picture of the experimental star simulator 



• All connections between the different components use 
FC/APC connectors to avoid parasitic back reflections. 

• The THT configuration has been optimized for the 
imaging of an unbalanced binary star with a high dynamics 
such as an exoplanet-star couple. Fo r this purpose we wil l 
use the theoretical results reported in I Armand et al.1 (|2008h 
which implies a redundant spacing of the array and an op- 
timized power distribution over the different telescopes to 
apodize the point spread function. 

• Our telescope array includes 8 apertures equivalent to 
a linear redundant configuration : the corresponding spatial 
frequencies sampling enables a convenient analysis of a com- 
plex linear object for a realistic experimental demonstration. 

The following sections summarize the THT bench structure. 
It consists of three main parts (cf Fig[5j): a star simulator, a 
telescope array and a combining interferometer. 



2.2 Technological options 

Our experimental set-up (see figure [5} has been designed 
and implemented thanks to the d ifferent skills devel- 
oped in our team for two decades (lAlleman et all f l995: 
ISimohamed fc Revnaudlll997l : iHuss et alll200ll : IPerrin et all 
120061 : lOlivier et all l2007h . Consequently, our THT experi- 
mental test bench uses optical fibres and couplers for the 
different optical functions to be implemented. However, we 
would like to stress that the use of guided optics compo- 
nents is not mandatory for the implementation of a THT. A 
classical design with classical components could be chosen 
if preferred. This point will remain a minor one as long as 
we will focus more on the demonstration of THT principle 
than on the technological aspects. The following items give 
the general framework of our experimental study. 

• The operating wavelength, all over this instrument, is 
A = 1, 55/xra to take advantage of the mature and available 
technologies of Telecom components. 

• Light propagation between the entrance pupil and the 
detector is achieved through monomode polarization main- 
taining fibres (panda fibres) . 



2.3 Star simulator 

The calibrated object is the first subsystem required for test- 
ing the imaging capability of a THT. For this first experi- 
mental demonstration, the selected astronomical target is a 
binary star with a convenient angular separation and ad- 
justable dynamics. 

For this purpose, the object consists of two tips of 
monomode Panda fibres glued on a V- groove. These 
monomode waveguides are fed by two independent Dis- 
tributed FeedBack lasers (DFB) with the same emitting 
wavelength and act as two incoherent point like sources. This 
way the object is spatially incoherent and the dynamics is 
controlled by adjusting the laser driving currents. A set of 
doublets and collimating lenses allows to provide an angular 
intensity distribution compatible with the spatial frequen- 
cies u sampled by our telescope array. In our experiment 
the angular separation #0, as seen by the telescope array, 
is 23.75/nrad. As our instrument is designed for a linear in- 
put polarization, a polarizing cube is inserted in the doublet 
spacing in order to select and fixes a linear vertical input 
polarization (not drawn on fig EJ) . The experimental setup 
can be seen in ng|6] 
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2.4 Design of the telescope array 

The telescope array arrangement has to be carefully selected 
to fit the sampling criter ia for a proper image analysis. As 
previously demonstrated (|Armand et all 20081 ) , high dynam- 
ics imaging capability requires a redundant array configu- 
ration. Consequently our telescope array must periodically 
sample the spatial frequency domain. The object dimension 
and the focal length of the collimator have to be determined 
by comparing the object spectrum and the spatial frequen- 
cies sampled by our instrument. The intensity I(At;y) ob- 
served in the image plane of the instrument is given by : 

I(At',y) = PSF®O(0) (23) 

where denotes the pseudo-convolution operator, PSF the 
point spread function of the instrument and O{0) , the object 
intensity distribution. In the Fourier domain, this relation- 
ship becomes a simple product: 

I{u) = T(u) x 6{u) (24) 

where 0(u) is the object intensity spectrum and T(u), the 
input pupil autocorrelation function. As shown on Fig. [Tithe 
periodicity of the spatial frequencies sampled by the tele- 
scope array has to be compliant with the classical Shannon 
sampling criterium. 

The smallest sampled frequency ( Br ™ n ) is related to 
the instrument field of view FV and the largest one (-^y 1 ) 
determines the instrument resolution R: 

FV = X/B mtn (25) 

R = \/B 

max (26) 

The telescope array resolution allows to discriminate the 
sharpest details to be observed on the object, and the field 
of view has to be adapted to the object's overall size. 
These design trade-offs lead to the following THT bench 



characteristics: 

FV — — 62firad 

J^min 



R — S. 9 fir ad 

-tJmax 

These characteristics have been chosen to be compatible 
with the observation of our laboratory binary star charac- 
terized by a 23, 75/j.rad angular separation. 

R < 6 < FV 

Notice that a redundant array, as proposed in this paper, in- 
duces a periodic PSF as reported on figure [8] The equation 
[25l gives the Field of View limits proper to avoid aliasing 
around the interferometric order. The area between two 
consecutive main lobes is called Clear Field (cf Fig(8}. In the 
case of polychromatic sources, only the zero interferometric 
order is achromatic and consequently appropriate for direct 
imaging. During the implementation of our breadboard we 
faced compactness constraint due to the mechanical dimen- 
sion of the launching assemblies. The focusing lenses are 
small enough to allow a compact linear array but the 3 axes 
nanopositioning mountings, required for the fine adjustment 
of the fibre inputs, cannot be aligned with a pitch smaller 
than 25 mm. So far, the real THT telescope array has been 
designed in a two-dimensional configuration but dedicated 
to observe a one dimensional object along the vertical axis. 
The projection of the telescope positions on the vertical axis 
is true to the expected linear array design as shown in Fig|9] 
In such a configuration, the instrumental response is exactly 
the same as with the linear telescope array when observing 
a vertical linear object. 

In order to image an unbalanced binary system with very 
high dynamics, the use of a set of suitable apodization co- 
efficients will optimize the PSF dynamics with a low decay 
of the resolution (cf Fig [13]). On each aperture, the control 
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Figure 7. Schematic analysis of the proper sampling criteria. The telescope array design has to be directly related to the angular 
dimension of the object: 4 ^ 
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of the intensity is achieved by means of mobile shutters (cf 
Fig JIOp actuated with a high position accuracy according t o 
the theoretical optimum distribution ([Armand et al.ll2008h . 
PSF dynamics D is denned as the ratio between the max- 
imum intensity on the PSF Imax and the higher intensity 
over the Clear Field imax (cf fiflU). 



Figure 9. Theoretical (white discs) and experimental (black 
discs) telescope array configuration: the projection of the tele- 
scope baseline along the vertical axis is true to the theoretical 
redundant configuration. These two configurations are equivalent 
to a ID vertical array when used with a ID vertical object. 



2.5 The interferometeric combiner 

The last part of the system is the optical field combiner mix- 
ing the contributions of the 8 telescopes (cf Fi £il2[) . Each in- 
terferometric ar m includes a fibre d elay line and an optical 
path modulator ([Olivier et afeOO/t ). The 8- fibre arms have 
been cut with a few mm accuracy in order to reduce the 
optical path differences as much as possible. The fibre de- 
lay lines allow to adjust the optical path with an accuracy 
of few \ivn. The fibre optical path modulators temporally 
reproduce the linear phase variation observed in the image 
plane of "classical" hypertelescopes (cf eq |13[) and allow the 
fine residual optical path compensation. A National Instru- 
ment virtual instrument and a voltage generator have been 



implemented and developed to drive the piezoelectric ac- 
tuators of the fibre optical path modulators. A set of high 
voltage amplifiers allows to reach a proper range of com- 
mand voltages. Experimentally, the optical path control is 
achieved with an optical path sensitivity over 0.01/ira. 
The optical fields emerging from the 8 interferometric arms 
reach a 8 to one polarization maintaining (PM) coupler to 
achieve the interferometric mixing. At the output an In- 
GaAs photodiode detects the interferometric signal that is 
recorded through a standard 12 bits ADC voltage acquisi- 
tion system. 
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Figure 10. THT input pupil of the telescope array. The lens 
doublets act as telescopes and mobile shutters allow to adjust the 
intensity ratio between the interferometric arms in order to get 
an apodized PSF. 
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Figure 11. PSF dynamics D is denned as the ratio between the 
maximum intensity on the PSF Imax and the higher intensity 
over the Clear Field imax 



3 EXPERIMENTAL TEST AND IMAGING 
CAPABILITIES OF THE THT TEST-BENCH 

3.1 PSF experimental results 

In order to calibrate the imaging properties of our instru- 
ment, the first investigation has been to characterize its 
point spread function. For this purpose, the telescope array 
has been illuminated by a plane wave using only a single 
point-like source ( i.e. switching on only 1 fibre of the ob- 
ject). A servo control system not available during this exper- 
iment, the 8 telescope cophasing has been achieved manually 
taking advantage of the relative stability of our instrument. 
The voltage offsets of the electronic commands, driving the 
optical path modulators, have been adjusted in order to in- 
crease the dynamics as well as possible. This is possible if 
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Figure 12. The THT interferometric beam combiner including 
the optical fibre delay lines and the optical path modulators. 
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Figure 13. Theoretical influence of the apodization process on 
the instrument PSF: the apodization allows to observe high dy- 
namics objects with a relatively low decay of the resolution (top 
linear scale, bottom log scale). 



you carrefully avoid any perturbation of the instrument over 
the cophasing process. For this purpose, the instrument has 
been acoustically and thermally baffled. The first experimen- 
tal results are shown in Fig |14l They are consistent with the 
theoretical PSF simulation proposed in Fig |131 The second 
main feature to be analyzed is the influence of the apodiza- 
tion process as shown on fig[T4l Using the ap odozation coeffi- 
cient computed in ref ([Armand et al ll2008l ) allows to reduce 
consequentially the ripple in the clear field all around the 



8 L. Bouyeron, S. Olivier and Other 



Experimental P&F without apodization ' 
Experimental PSF with apodization 




-30 



-20 -10 10 20 
Angular position (^rad) 



30 



10' 



10° 



>* 10" 



10"' 



10"~ 



Experimental PSF without apodization 
Experimental PSF with apodization 





1.0E+3 1.0E+4 

rejection rate 



Figure 16. Histogram on minimum dynamics as function of fluc- 
tuations on the apodization coefficient: this histogram is obtained 
from one million of simulated PSF. For each PSF, the apodization 
coefficient values are determined randomly using the experimen- 
tal statistic ( mean value and RMS). This numerical simulation 
demonstrates that in the worst case of 1500 dynamics is far better 
than the experimental result. Accordingly, this error source has a 
minor effect. 
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Figure 14. Experimental PSF obtained on our 8-telescope THT 
test bench: The apodization process allows to reach a 300 dy- 
namics between the maximum and the clear field area, (top linear 
scale, bottom log scale) 
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Figure 15. PSF dynamics as function of time: example of the 
evolution of experimental PSF dynamics over 180 s. 



PSF peak. The best dynamics experimentally observed with 
our breadboard are in the range of 300 but as shown on fig- 
ure [15] these results are not stable over time, probably due to 
environmental perturbations such as acoustic waves, vibra- 
tions or temperature fluctuations. Nevertheless, such results 
have been obtained reliably and give a first indication on the 
THT imaging potential capabilities. As we will demonstrate 
in the next paragraph, the major limitation on the dynam- 
ics results from the manual cophasing that remains of poor 
quality when compared with a servo controlled system. The 



main effect is a reduction of the dynamics that is particu- 
larly clear in the apodized configuration. These first tests 
demonstrate the possibility to image a point-like source on 
axis. 

3.2 Analysis of experimental results 

The phenomena able to reduce the PSF dynamics can be 
listed as follows: 

• The residual optical path difference between the 8 in- 
ter ferometric arms. 

• The differential polarization properties between the 8 
fields to be mixed due to differential birefringent properties 
of the interferometric arms. 

• The difference betw een the theoretical results reported 
in lArmand et all (|2008h and the apodization really applied 
in the experiment. 

• The differential dispersion properties between the 8 in- 
terferometric arms. 

This last point is not relevant in our experiment due to the 
use of a quasi-monochromatic source. 

The influence of the third point has been analyzed by mod- 
eling the effect of the experimental uncertainties on apodiza- 
tion coefficients. According to the error experimentally ob- 
served, the en coefficients are randomly changed in a nu- 
merical simulation and the final PSF is computed. Table [1] 
summarizes the experimental uncertainties observed in our 
test bench. This statistical distribution on ai coefficient has 
been used as input for our numerical simulations. A set of 
PSF has been computed varying the ai allowing to plot a 
histogram of the rejection ratio ( Fig. [16]). This curve demon- 
strates that the probability to get dynamics lower than 1500 
is practically zero. This dynamics is far higher than the ob- 
served value. Accordingly, the experimental uncertainty on 
the apodization coefficient cannot explain the observed ex- 
perimental rejection ratio and cannot be identified as re- 
sponsible for the 300 experimental dynamics. 
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Arm 


apodization 


theoretical 


experimental 


i 


coeff 


normalized 


intensity 




CLi 


intensity 


range 


1 - 8 


0.104 


0.011 


[0.005;0.016] 


2 - 7 


0.364 


0.133 


[0.122;0.143] 


3 - 6 


0.724 


0.525 


[0.496;0.553] 


4 - 5 


1 


1.000 


[0.985;1.015] 



Table 1. Apodization parameters 



spatial 
frequencies 


experimental 
contrast 


1 


0.9996 (±4.10- 4 ) 


2 


0.9995 (±4.10- 4 ) 


3 


0.9993 (±3.10- 4 ) 


4 


0.9994 (±4.10- 4 ) 


5 


0.9992 (±3.10- 4 ) 


6 


0.9992 (±3.10- 4 ) 


7 


0.9991 (±3.10- 4 ) 



Table 2. Experimental contrasts as function of the sampled spa- 
tial frequencies achieved with pairs of telescope. 



In order to carry on our investigation on with the exper- 
imental origin of the limited dynamics, the next point under 
study is the polarization behaviour. The optical fibre arms 
of the interferometer are made of polarization-maintaining 
monomode fibres. This choice results from the necessity to 
avoid any loss of the polarization coherence between the op- 
tical fields to be interferometrically mixed. There can be two 
problems : 

• Input linear polarization can be misaligned with the 
used neutral axis at the input of waveguides. 

• Crosstalk between the two polarization modes corre- 
sponding to the two neutral axes. 

In both cases, it results in an incoherent background in- 
tensity propagating along the undesirable polarization axis. 
This parasitic intensity leads to a background limit reducing 
the dynamics to be reached by the instrument. In our exper- 
iment, care has been taken to limit this problem by means 
of polarizers placed at the entrance pupil and after the 8 to 
1 optical coupler. Our bread board has been characterized 
by measuring the contrast with each telescope couple. All 
the contrasts remain better than 0.999 as reported in Ta- 
ble|2] According to these experimental data, it is possible 
to compute a global PSF introducing a random polariza- 
tion crosstalk for each interferometric arm. It results in a 
histogram of the rejection ratio to be reached according to 
such a perturbation fFig J17[) . This simulation shows that the 
probability to get a rejection ratio lower than 500 is prac- 
tically zero. Accordingly, the polarization behaviour cannot 
be identified as responsible for the 300 dynamics observed 
in our experiment. Consequently, the 300 dynamics can only 
be explained by the residual optical path differential fluctu- 
ations between the 8 interferometric arms. 



3.3 THT imaging capabilities 

The second stage of our experimental investigation has been 
to check the invariance of this point spread function over the 
imaging field of view. As a preliminary test, we moved the 




1,0E+3 1.0E+4 

rejection rate 



Figure 17. Birefringence rejection ratio limitation: This his- 
togram obtain by 100000 PSF simulation shows the birefringence 
effect on the rejection ratio in the case of our experimental setup 
(contrast = 0.999). In the worst case, the rejection ratio is limited 
to 500. 
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Figure 18. PSF translation quasi invariance: Moving the point- 
like source in THT field of view allows to observe global shift of 
the PSF. The repeatability is experimentally demonstrated only 
at a 100:1 level due to the turbulence phase fluctuation. 



point like source in the object plane by shifting the optical 
fibre tip using manual translations. As shown on figure fTHl 
the PSF was simply shifted with a constant width of the 
main lobe. The intensity on the clear field shows some fluc- 
tuations leading to dynamics losses. These imperfections are 
due to environmental parameter fluctuations resulting from 
the lack of a real time cophasing system. This demonstrates 
that our instrument is linear and translation invariant and 
therefore suitable for imaging . 
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Figure 19. Experimental image of an unbalanced binary system 
with a 20 power ratio using our 8-telescope array in the apodized 
configuration(top linear scale, bottom log scale). 



In order to illustrate the imaging capabilities of our instru- 
ment, we have used a binary star as source. The two fibre 
tips of the laboratory object receive the light from two dif- 
ferent lasers operating at the same wavelength. Adjusting 
the driving current allows to simulate any unbalanced bi- 
nary source. In the absence of a real time cophasing system, 
we used a secondary star intensity higher than the clear field 
intensity fluctuations. 

The experimental result is shown in fig|19land exhibits 
a main lobe surrounded by two side lobes. The bright one 
is related to the main star and the two others to its com- 
panion with the interference order on the right and the 
first one on the left of this figure. Notice that in the case of 
polychromatic sources, only the zero interferometric order 
will be used as being achromatic. The source characteristics 
with a 23. 75 fir ad separation and a ~ 20 intensity ratio are 
correctly recovered in the image. Results are consistent with 
the theoretical convolution between the point spread func- 
tion and the object intensity distribution according to eq|23l 
Consequently, over the instrument clear field the instrument 
imaging capabilities have been demonstrated. 



3.4 Preliminary 2D PSF observation 

This first experimental study was mainly dedicated to the 
demonstration of the imaging capabilities of a temporal hy- 
pertelescope on a one dimensional object. Due to mechanical 
constraints, the telescope array has been designed in a 2D 
configuration. Even if this array is not absolutely optimized 
for a two dimensional imaging, it is possible to test the 2D 
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Figure 20. Preliminary test on the 2D imaging capabilities of a 
Temporal Hypertelescope. Comparison between ideal (top) and 
experimental (middle) 2D PSF obtained with our THT bench 
(bottom)absolute value of the difference between the ideal PSF 
and the experimental image. 



PSF to demonstrate the operation of the 2D imaging pro- 
cess. This picture can be obtained row by row like in a TV 
raster scan using the phase modulations reported in equa- 
tion eqJTl A set of phase shifts ipm is sent through the driving 
voltage electronics taking into account the projection of the 
telescope i baseline along the y direction : 



<Pio(yo) = 27T.Pi.yo 



(27) 



This equation describes the phase shift that could be ob- 
served for telescope i in the image plane along the verti- 
cal axis in a spatial configuration. As demonstrated in ref 
(jRevnaud fe Delagd (j2007h ) it is possible this way to get a 
2D information using a temporal hypertelescope. We tested 
this process by recording the 2D PSF using our experimen- 
tal test bench. The theoretical and experimental patterns 
are reported on Fig |20l The array cophasing was manually 
achieved during the first scan on the central row (yo = 0). 
After this initialization process, it is possible to scan all the 
field of view taking advantage of the relative stability of the 
telescope array. This full operation lasts about 20 seconds to 
get 50 rows over the 2D PSF. In order to quantify the dis- 
crepancies between ideal (Figure [20] top) and experimental 
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(Figure [20] middle) 2D PSF, Figure [20] bottom reports the This paper has been typeset from a TgX/ I^I^X file prepared 
absolute value of the difference between the two images. The by the author, 
significant fluctuations are probably due to the lack of an ac- 
tive phase control system during the 2D scan. A fine study 
of this limitation will be achieved after the implementation 
of a servo control system of our test bench. Nevertheless, 
these results are promising and could be enhanced with a 
convenient cophasing system. 



4 CONCLUSION 

The aim of this study was to experimentally demonstrate 
the operation of a temporal hypertelescope. This experi- 
mental demonstration was achieved by testing the instru- 
ment Point Spread Function as a first stage. The second 
experiment demonstrated the possibility to observe an un- 
balanced binary star with a 20 flux ratio. During this first 
test, the PSF dynamics was limited to 300 due to the lack 
of servo control system. The next stage of this study will be 
the design and implementation of a servo control system to 
accurately cophase our telescope array in order to enhance 
the dynamics of our instrument. 
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